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ABSTRACT 

The structure of e + e~ — ► bbg events was studied using Z° decays recorded in the 
SLD experiment at SLAC. Three-jet final states were selected and the CCD-based 
vertex detector was used to identify two of the jets as b or b. Distributions of the 
gluon energy and polar angle were measured over the full kinematic range for the first 
time, and compared with perturbative QCD predictions. The energy distribution is 
potentially sensitive to an anomalous b chromomagnetic moment k. We measured k 
to be consistent with zero and set the first limits on its value, —0.17 < k < 0.11 at 
95% c.l. 
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The observation of e + e~ annihilation into final states containing three hadronic jets, 
and their interpretation in terms of the process e + e~ — > qqg Q, provided the first 
direct evidence for the existence of the gluon, the gauge boson of the theory of strong 
interactions, Quantum Chromodynamics (QCD). In subsequent studies the jets were 
usually energy ordered, and the lowest-energy jet was assigned as the gluon; this is 
correct roughly 80% of the time, but preferentially selects low-energy gluons. If the 
gluon jet could be tagged explicitly, event-by-event, the full kinematic range of gluon 
energies could be explored, and more detailed tests of QCD could be performed M. 
Due to advances in vertex-detection this is now possible using e + e~ — ► bbg events. The 
large mass and relatively long lifetime, ~ 1.5 ps, of the leading B hadron in 6-quark 
jets || lead to decay signatures which distinguish them from lighter-quark (u, d, s or 
c) and gluon jets. We used our 120-million-pixel CCD vertex detector || to identify 
in each event the two jets that contain the B hadrons, and hence to tag the gluon jet. 
This allowed us to make the first measurement of the gluon energy and polar-angle 
distributions over the full kinematic range. 

Additional motivation to study the bbg system has been provided by measurements 
involving inclusive Z° — > bb decays. Several reported determinations || of Rt, = 
T(Z° — >bb)/T(Z° —>qq) and the Z°-b parity-violating coupling parameter, Aj,, differed 
from Standard Model (SM) expectations at the few standard deviation level. Since 
one expects new high-mass- scale dynamics to couple to the massive third-generation 
fermions, these measurements aroused considerable interest and speculation. We have 
therefore investigated in detail the strong-interaction dynamics of the 6-quark. We 
have compared the strong coupling of the gluon to 6-quarks with that to light- and 
charm-quarks 0, as well as tested parity (P) and charge©parity (CP) conservation at 
the bbg vertex [[?]. Here we study the structure of bbg events via the distributions of the 
gluon energy and polar angle with respect to (w.r.t.) the beamline. We compare these 
results with perturbative QCD predictions, including a recent calculation at next-to- 
leading order (NLO) which takes quark mass effects into account M. 
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In QCD the chromomagnetic moment of the b quark is induced at the one-loop level 
and is of order a s /7T. A more general bbg Lagrangian term with a modified coupling |§ 
may be written: 

C hl9 = gsbT a {^ + —^——(k - iR lh )}bG», (1) 

where k and k parameterize the anomalous chromomagnetic and chromoelectric mo- 
ments, respectively, which might arise from physics beyond the SM. The effects of the 
chromoelectric moment are sub-leading w.r.t. those of the chromomagentic moment, 
so for convenience we set k to zero. A non-zero k would modify || the gluon energy 
distribution in bbg events relative to the standard QCD case. Hence we have used our 
data to set the first limits on k. 

We used hadronic decays of Z° bosons produced by e + e~ annihilations at the SLAC 
Linear Collider (SLC) which were recorded in the SLC Large Detector (SLD) |TU|. The 
criteria for selecting Z° decays, and the charged tracks used for flavor-tagging, are 
described in || |TH| . We applied the JADE algorithm |L2| to define jets, using a scaled- 
invariant-mass criterion y cut = 0.02. Events classified as 3-jet states were retained if 
all three jets were well contained within the barrel tracking system, with polar angle 
| cos8j et \ < 0.71. From our 1993-95 data samples, comprising roughly 150,000 hadronic 
Z° decays, 33,805 events were selected. In order to improve the energy resolution the 
jet energies were rescaled kinematically according to the angles between the jet axes, 
assuming energy and momentum conservation and massless kinematics [0]. The jets 
were then labelled in order of energy such that Ei > E 2 > E 3 . 

Charged tracks with a large transverse impact parameter w.r.t. the measured in- 
teraction point (IP) were used to tag bbg events 0. The resolution on the impact pa- 



rameter, projected in the plane normal to the beamline, d, is cr^ = 11©70 / \p±Vsm 6) 
/im, where p± is the track transverse momentum in GeV/c, and 9 the polar angle, 
w.r.t. the beamline. The flavour tag was based on the number of tracks per jet, N 3 S , 



et 
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with d/aa > 3. Events were retained in which exactly two jets were 6-tagged by re- 
quiring each to have > 2, and in which the remaining jet had N 3 J? < 2 and was 
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hence tagged as the gluon; 1329 events were selected. The efficiency for selecting true 
bbg events is 8.3%. This was estimated using a simulated event sample generated with 
JETSET 7.4 [|n|, with parameter values tuned to hadronic e + e~ annihilation data [14 



combined with a simulation of S-decays tuned to T(4S) data IIJ and a simulation of 
the detector. The efficiency peaks at about 11% for 15 GeV gluons. Lower-energy 
gluon jets are sometimes merged with the parent 6-jet by the jet-finder. At higher 
gluon energies the the correspondingly lower-energy 6-jets are harder to tag, and there 
is also a higher probability of losing a jet outside the detector acceptance. 

For the selected event sample, Fig. 1 shows the N J Jg distributions separately for jets 
1, 2 and 3. In about 15% of cases the gluon-tagged jet is not the lowest-energy jet (jet 3). 
The simulated contributions from true gluons are indicated, and the estimated gluon 
purities [[Uj are listed in Table |l| The inclusive gluon purity of the tagged-jet sample 
is 95%. With this sample we formed the distributions of two gluon-jet observables, 
the scaled energy x g = 2E g \ non / y/s, and the polar angle w.r.t. the beamline, 9 9 . The 
distributions are shown in Fig. 2. The simulation is also shown; it reproduces the data. 

The backgrounds were estimated using the simulation and are of three types: non-66 
events, 66 but non-bbg events, and true bbg events in which the gluon jet was mis-tagged 
as a 6-jet. These are shown in Fig. 2. The non-66 events (~ 5% of the bbg sample) 
are mainly ccg events, 90% of which had the gluon correctly tagged. There is a small 
contribution (~ 0.1% of the bbg sample) from light-quark events. The dominant back- 
ground is formed by 66 but non-bbg events. These are true 66 events which were not 
classified as 3-jet events at the parton level, but which were mis-reconstructed and 
tagged as 3-jet bbg events in the detector using the same jet algorithm and y cut value. 
This arises from the broadening of the particle flow around the original 6 and 6 di- 
rections due to hadronisation and the high-transverse-momentum S-decay products, 
causing the jet-finder to reconstruct a 'fake' third jet, which is almost always assigned 
as the gluon. The population of such fake gluon jets peaks at low energy (Fig. 2(a)), 
as expected. Mis-tagged events comprise less than 1% of the bbg sample. 
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The distributions were corrected to obtain the true gluon distributions D true (X) by 
applying a bin-by-bin procedure: D true (X) = C(X) (D raw (X) - B(X)), where X = x g 
or cos6 g , D raw (X) is the raw distribution, B(X) is the background contribution, and 
C(X) = D^ u c (X) / D M Q n (X) is a correction that accounts for the efficiency for accept- 
ing true bbg events into the tagged sample, as well as for bin-to-bin migrations caused 
by hadronisation, the resolution of the detector, and bias of the jet-tagging technique. 
Here -D^c(X) is the true distribution for MC-generated bbg events, and D M Q n (X) is 
the resulting distribution after full simulation of the detector and application of the 
same analysis procedure as applied to the data. 

As a cross-check, an alternative correction procedure was employed in which bin- 
to-bin migrations, which can be as large as 20%, were explicitly taken into account: 
D true (Xi) = M{X h Xj) {D raw {X j )-B{X j ))/e{X i ), with the unfolding matrix M(X U X 5 ) 
defined by .D^f(Xj) = M(Xi, Xj)D M c n (Xj), where true bbg events generated in bin 
i may, after reconstruction, be accepted into the tagged sample in bin j. e(X) is 
the efficiency for accepting bbg events in bin i into the tagged sample. The result- 
ing distributions of x g and cos^ are statistically indistinguishable from the respective 
distributions yielded by the bin-by-bin method. 

The fully-corrected distributions are shown in Fig. 3. Since, in an earlier study ||, 
we verified that the overall rate of bbg-event production is consistent with QCD expec- 
tations, we normalised the gluon distributions to unit area and we study further the 
distribution shapes. The x g distribution rises, peaks around x g ~ 0.15, and decreases 
towards zero as x 9 — > 1. The peak is a kinematic artifact of the jet algorithm, which 
ensures that gluon jets are reconstructed with a non-zero energy which depends on the 
y c value. The cos^ distribution is flat. 

We have considered sources of systematic uncertainty that potentially affect our 
results. These may be divided into uncertainties in modelling the detector and un- 
certainties in the underlying physics modelling. To estimate the first case we sys- 
tematically varied the track and event selection requirements, as well as the tracking 
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efficiency || [TT[ . In the second case parameters used in our simulation, relating mainly 
to the production and decay of charm and bottom hadrons, were varied within their 
measurement errors [[UJ. For each variation the data were recorrected to derive new x g 
and cos8 g distributions, and the deviation w.r.t. the standard case was assigned as a 
systematic uncertainty. None of the variations affects our conclusions. All uncertain- 
ties were conservatively assumed to be uncorrelated and were added in quadrature in 
each bin of x g and cos9 g . 

We compared the data with perturbative QCD predictions for the same jet al- 
gorithm and y c value. We used leading-order (LO) and NLO results based on recent 
calculations || in which quark mass effects were explicitly taken into account; a 6-mass 
value of m b {m z ) = 3GeV/c 2 was used [[T7[]. We also derived these distributions using 
the 'parton shower' (PS) implemented in JETSET. This is equivalent to a calculation 
in which all leading, and a subset of next-to-leading, lny c terms are resummed to all 
orders in a s . In physical terms this allows events to be generated with multiple orders 
of parton radiation, in contrast to the maximum number of 3 (4) partons allowed in 
the LO (NLO) calculations, respectively. Configurations with > 3 partons are relevant 
to the observables considered here since they may be resolved as 3-jet events by the 
jet-finding algorithm. 

These predictions are shown in Fig. 3. The calculations reproduce the measured 
cos# 9 distribution, which is clearly insensitive to the details of higher-order soft parton 
emission. For x g , although the LO calculation reproduces the main features of the 
shape of the distribution, it yields too few events in the region 0.2 < x g < 0.5, and too 
many events for x g < 0.1 and x g > 0.5. The NLO calculation is noticeably better, but 
also shows a deficit for 0.2 < x g < 0.4. The PS calculation describes the data across 
the full x g range. The \ 2 f° r the comparison of each calculation with the data is given 
in Table ^. These results suggest that multiple orders of parton radiation need to be 
included, in agreement with our earlier measurements of jet energy distributions using 
flavor- inclusive Z° decays |]18[ . We also investigated LO and NLO predictions based on 
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matrix elements implemented in JETSET which assume massless quarks. The resulting 
distributions are practically indistinguishable from the massive ones, even though the 
large 6-mass has been seen ]17] to affect the bbg event rate at the level of 5%. The 
effect of varying a s within the world-average range is similarly small. 

We conclude that perturbative QCD in the PS approximation accurately repro- 
duces the gluon distributions in bbg events. However, it is interesting to consider the 
extent to which anomalous chromomagnetic contributions are allowed. The Lagrangian 
represented by Eq. |I] yields a model that is non-renormalisable. Nevertheless tree-level 
predictions can be derived and used for a 'straw man' comparison with QCD. For 
illustration, the effect of a large anomalous moment, k = 0.75, on the shape of the 
x g distribution is shown in Fig. 3(a); there is a clear depletion of events in the re- 
gion x g < 0.5 and a corresponding enhancement for x g > 0.5. By contrast the shape 
of the cos6 g distribution is relatively unchanged (not shown), even by such a large k 
value. In each bin of the x g distribution, we parametrised the leading-order effect of 
an anomalous chromomagnetic moment and added it to the PS calculation to arrive 
at an effective QCD prediction including the anomalous moment at leading-order. A 
X 2 minimisation fit was performed to the data with k as a free parameter, yielding 
k = —0.029 ± 0.070(stat.)^QQQ3(syst.), which is consistent with zero within the errors, 
with a x 2 °f 9-3 for 9 degrees of freedom. The distribution corresponding to this fit 
is indistinguishable from the PS prediction (Fig. 3(a)) and is not shown. Our result 
corresponds to 95% confidence-level (c.l.) upper limits of —0.17 < k < 0.11. 

In conclusion, we used the precise SLD tracking system to tag the gluon in 3- 
jet e + e~ — > Z° — > bbg events. We studied the structure of these events in terms of 
the scaled gluon energy and polar angle, measured for the first time across the full 
kinematic range. We compared our data with perturbative QCD predictions, and 
found that the effect of the 6-mass on the shapes of the distributions is small, that 
beyond-LO QCD contributions are needed to describe the energy distribution, and 
that the parton shower prediction agrees best with the data. We also investigated an 
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anomalous 6-quark chromomagnetic moment, k, which would affect the shape of the 
energy distribution. We set 95% c.l. limits of —0.17 < k < 0.11. As far as we are 
aware, these are the first such limits on an anomalous quark chromomagnetic coupling. 
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Jet label 


# Tagged gluon jets 


Purity 


3 


1140 


96.1 % 


2 


155 


90.2 % 


1 


34 


65.7 % 



Table 1: Estimated purities of the tagged gluon-jet samples. 



QCD Calculation 


X 2 : a; s (10 bins) 


LO m b (m z ) = 3 GeV/c 2 


73.6 


NLO m b (m z ) = 3 GeV/c 2 


24.3 


PS M 6 = 5 GeV/c 2 


9.5 



Table 2: \ 2 f° r t ne comparison of the QCD predictions with the corrected data. 
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Figure 1: The iV^ distributions for jets in 66g-tagged events, labelled according to jet 
energy (dots); errors are statistical. Histograms: simulated distributions showing jet 
flavour contributions. 
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Figure 2: Raw measured distributions of (a) x g and (b) cos# 9 (dots); errors are statis- 
tical. Histograms: simulated distributions including background contributions. 
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Figure 3: Corrected distributions of (a) x g and (b) cos# 9 (dots); errors are statistical. 
Perturbative QCD predictions (see text) are shown as lines joining entries plotted at 
the respective bin centers. 
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